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Abstract 

We discuss a reaction-diffusion model in one dimension subjected to an external driving force. Each 
lattice site may be occupied by at most one particle. The particles hop with rates (1 ± »7)/2 to the right 
or left nearest neighbour site if it is vacant, and annihilate with rate one if it is occupied. We compute 
the long time behaviour of the space dependent average density in states where the initial density profiles 
are step functions. We also compute the exact time dependence of the particle density for uncorrelated 
random initial conditions. The representation of the uncorrelated random initial state and also of the 
step function profile in terms of free fcrmions allows for the calculation of time- dependent higher order 
correlation functions. We outline the procedure using a field theoretic approach. Finally, we show how 
this gives rise to predictions on experiments in TMMC exciton dynamics. 

PACS numbers: 05.40.-f-j, 05.50.-fq, 02.50.Ga, 82.20.Mj 
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1 Introduction 

In recent years there has been much interest in one-dimensional reaction-diffusion systems. Despite their 
simpUcity, they exhibit a very rich dynamical behaviour and have been used as models for experimental 
systems [0. Among them, the study of diffusion of hard-core particles (e.g. with exclusion) plays a central 
role j|, H, ^. Such systems are usually described in terms of classical (mean field) equations. However, in low 
dimensional systems, where the diffusion of the reactants is not sufficiently rapid compared to the reaction 
rate, large fluctuations persist and lead to a breakdown of the mean field approximation. Thus more rigorous 
methods are needed, e.g. a description of the model in terms of a master equation. 

In this paper we will consider a model which allows for the full solution of the master equation. It is 
defined on a ring of L sites where each lattice site may be occupied by at most one particle. This exclusion 
principle models a hard-core repulsion of the particles. These particles (which will be denoted by A) hop 
with rates (1 ±7y)/2 to the right or left nearest neighbouring sites respectively, if the site is vacant (denoted 
by 0), and annihilate with rate A if it is occupied: 

A0 — > 0yl (l + 77)/2 (1) 

%A ^ A% (l-?7)/2 (2) 

A A — > 00 A. (3) 

We will study A = 1. This case is the one which is fully solvable, and besides the results obtained are of 
experimental relevance for the study of exciton dynamics in TMMC polymers (see below). This model, which 
is closely related to zero temperature Glauber dynamics has already been studied by many authors, see 
e.g. i - 0. 

It turns out to be convenient to write the master equation of the process in terms of a quantum spin 
Hamiltonian. Since the particle model involves two states per site, it maps into a spin 1/2 system, the 
evolution operator of which is the quantum spin Hamiltonian. The Jordan- Wigner transformation allows 
this Hamiltonian to be written in terms of fermion operators and reduces for A = 1 to a free fermion system. 
This somewhat surprising simplification allows for the exact calculation of many dynamical properties of the 
system. In particular it has been known for a long time that the density for an initially completely filled 
state decays as t^^/"^ |6| corresponding to a diffusive "critical" behaviour. 
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We extend these results and calculate explicitly the time evolution of the density for uncorrelated random 
initial states with any given density < p < 1 . We also show how to calculate higher order correlation func- 
tions. Moreover, we study the time evolution of the space dependent average density for non-translationally 
invariant initial states, namely the step function and half step function in the density profile. With this we 
obtain the development in time of the reaction edge of a system which is initially completely empty to the 
left of the origin and full (half-full) to the right of it. Also, we have found a convenient way to calculate 
time-dependent higher order correlation functions for the step function initial state. 

The interest of studying this problem is not merely theoretical. This model can be mapped by a similarity 
transformation into the coagulation problem where we have the reaction {A A — > A , $ A) rather than 
the pair annihilation A A — > It has been suggested |l^ to use this lattice model for the 

description of exciton kinetics in the carrier suhsUnce N{CH3)iMnCl3 (TMMC) §: Excitons of the Mn2+ 
ion are generated after pulsed laser excitations of the sample and diffuse along the widely separated MnCl^ 
chains. Experimental data show that a single exciton has a decay time of about 0.7ms. The on-chain 
hopping rate is 10^^ — lO^^s^^. If two excitons arrive on the same Mn?^ ion, they undergo a coagulation 
reaction {A-\- A — > A), under the emission of radiation. The reaction time for this process is assumed to 
be ~ lO^^'^s which is much faster than the diffusion time. Measurements of the time resolved luminescence 
in TMMC, in a time scale of a picosecond to a millisecond allow the observation of the decay of the exciton 
population. Two characteristics of this system are especially relevant to the lattice coagulation model we 
are using to describe it. Firstly, given that the reaction process is much faster than the diffusion process, 
the coagulation rate for the lattice model is determined to be equal to the hopping rate: If two particles (i.e. 
excitons) are on neighbouring lattice sites (i.e. Mn^+ ions), either of them will hop onto the site occupied by 
the other with the (unit) diffusion rate. Almost instantaneously one of them will annihilate and as a result 
one particle will be left, with equal probability on either site. Secondly, the single exciton decay is very slow 
compared with the other processes and we can, as a first approximation to the dynamics, neglect it. The 
model then exactly maps into the annihilation problem with A = 1. Hence, we can obtain expressions for 
physical quantities that can be compared with experiments. 

The structure of this paper is as follows: In section 2 we introduce the description of this system in 
terms of the master equation and analyse its mapping to the quantum spin problem. In section 3, we sketch 
how this problem can be converted into a problem of spinless fermions, by means of the Jordan- Wigner 
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transformation, and how to calculate the relevant quantities in this fermion language. We proceed to the 
calculation of the density for the uncorrelated states, in section 4. We postpone to Appendix A the techniques 
which allow the calculation of higher order correlation functions. In section 5 we present the calculation of 
the average space dependent density for the step and half step function states, and again leave the details 
related to the calculation of higher order correlation functions to an appendix, in this case Appendix B. 
We repeat the former calculation for a half step function state. In section 6, we briefly review the concept 
of similarity transformation, and use it to calculate quantities which can be compared with experiments. 
Finally, in section 7, we present our conclusions and some open questions. 



2 The master equation and the mapping to the spin system 

Our system is defined on a lattice of L sites with periodic boundary conditions. Each site has assigned 
a particle occupation number rij = 0, 1 (1 < J < L). The set n — {ni, . . . ,7^L} characterizes a given 
configuration of particles. The process is completely defined by the probability distribution P{n, t) which 
gives the probability for the system to be in the configuration n at time t. P{n, t) satisfies a master equation 
which is conveniently expressed in terms of a quantum Hamiltonian language . 

The idea is to assign to each of the 2^ configurations of this system a normalized basis vector | n ) in a 
Hilbert space. The probability distribution is then mapped to a state vector 

\P{t)) = E nn,t)\n). (4) 

Using this notation it is then possible to write the master equation in a "Schrodinger form" 

||P(t)) = -H\P{t)) (5) 

where _ff is a linear, in general non-hermitian operator, containing all the information on the stochastic 
process. Equation (j^) has the formal solution 

\P{t)) = exp(-i/t) |P(0)). (6) 

Now one can express the time dependent averages of any function A(ni, . . . , n^) of the occupation num- 
bers as 

= {s\A{n^,...,nL)\P{t)) (7) 



where the n^ 's are the projection operators on states with a particle on site j of the chain and ( s | is defined 



as 



(8) 



nex 



Using (H) one sees that 



s\H = 0. 



(9) 



This expresses the conservation of the total probability. 

Since we are working with a two state model the projection operators hj can be represented in terms of 
Pauli matrices 



1 



'00' 



yO Ij 



(10) 



Also, for convenience, we introduce the operators sf — {a^ ± ia^ )/2. In the convention used just above 



(11) 



creates a particle at site j, and 



/o ^ 



(12) 



annihilates a particle at site j. 

With these operators the quantum Hamiltonian which describes the stochastic dynamics of the pair 
annihilation process is given by 



i 1 



+ 2 (1 ~ ^) [ (1 - "j+i - sj ; 



+ X{hjhj+i - s+J+^J}. 



(13) 



We conclude this section noting that we can change from the Schrodinger representation of the operators 
to the Heisenberg one in a similar way as in quantum mechanics. The time dependent operators are defined 

by 



A{t) = exp(iJi) A cxp (-Ht) 



(14) 



and obey the equation of motion 
dA{t) 



dt 



H,A{t)]. (15) 



3 The fermion representation of the operators 

We introduce a fermion representation of the Hamihonian by means of the Jordan- Wigner transformation 
[|9|. Defining 



Q, = n (16) 

i=l 

a] = S- Q,-i (17) 
a, = (18) 
it is a simple exercise to show that the set of a| 's, a/ 's obey the fermion anticommutation relations 
{cLj , di } = { a] , a| } = 

{a,, a|} = S,,i. (19) 

Note that because of the periodic boundary conditions for the spin operators one has a]^^i — Ql 
and aL+i = Ql^i- Ql may be written Ql = (—1)^, where N — fij is the number operator. Since 
the action of H changes the particle number only by units of two, Qi, is a constant of motion and splits 
the Hilbert space into a sector with an even number of particles {Ql = 1) and into a sector with an odd 
number of particles {Ql = —1)- The application of aL+i into states belonging to one of this sectors gives 
ciL+i \ X) — — ai \ X) \i X e even sector and ol+i | -'^) — di \ X) \i X e odd sector. This implies that 
we have antiperiodic (periodic) boundary conditions for the even (odd) sector. Applying the equations ( p^ 
- ( |T^ ) to the Hamiltonian (|l3|), one obtains [Q 

L 
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F = - i ^{d]j^^aj + d]aj+i + 2aj+idj - 2 ) - ( A - 1 ) ^(aj+iflj - fijUj+i) 
i=i i=i 

L 

{hj = ajflj). Notice that when we take A = 1, the two body term vanishes, making this Hamiltonian a 
free fermion one. 
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We now define the Fourier transforms 

L 



Up 



hi 



e 4 

7i 



e 4 



a. 



(21) 



(22) 



where p runs over aU integers p = — -^j . . . , 0, . . . , — 1 in the sector with odd number of particles and 



over the half-odd integers p = —h — \, 



1 3 
' 2 ' 2 ' ■ 



^ — i in the even sector (we take L even). Also, the 



transformation ( pl| ) - ( p2[ ) is canonical and one has {bp.bq} — { , } — and { , &J } = <5p,q. In 



terms of the 6p's, foj's, the Hamiltonian (20) reads |T4| 



2tip \ 



. f '2ttp\ 
sm 



1 - cos ( -j^ I I 6l,6p + sin [~J^] b-pbp 



i'^E 



sin ^p^P- 



(23) 



Since is quadratic in the 5p's, fej's, the equations of motion (15) for these operators are linear and 
become easily solvable Mm. Their solution is 



bp{t) 



cot 



(1 _ e-(^''+--)*)L 



(24) 
(25) 



with 



= 1 — cos 



/ 27rp \ . . / 2t:p 

— 1 77 sm 

\ L J ' \ L 



(26) 



Having obtained the dynamical solution of the problem, we now turn to the description of the initial 
conditions in this fcrmion language. With the convention we adopted in section 2, the state with all spins up 
is the empty state, i.e. the vacuum. A state with particles in sites ji < i2 < ■ • ■ < Jn {N < L) is represented 

by 



Jl,j2, ■■■,JN , 



(27) 



As ji < j2 < ■ ■ ■ < jN, Qji-1 can be commuted with all the remaining Sj operators and applied to the 
vacuum, giving one. The same holds for the remaining Q operators, giving 



I ji, j2, ■■■Jn) = a], a]^ ... a]^ I 0) 



(28) 
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i.e. the states are just Slater determinants in the position representation. Also, from the condition 
{s \H = 0, one obtains 

{s\ = (0 I + cot (^) 6pLp) + (0 |6o + cot ('^') V6-P') 

= ( I exp cot ( ^ ) bpb^p) + {0 \ bo exp *^°M ) ^p' ^~p' ^ 
p p' ^ ' 

= (s I™™ + (s r'*'*. (29) 

The product and sum over p, p' run over p — — ^ (even sector) and over p' = 0, . . . , — 1 (odd 

sector), respectively and [s \ obey (s | ^ ' ii = separately. 

For simplicity, we only consider states belonging to the even sector of the Hilbert space. Since the 
calculation of all correlation functions implies the action of an even number of fermion operators to such a 
state, and then the contraction with (s |, one can replace (s | just by (s ^'^""^ because the contraction with 
( s 1°'*'* gives 0. One property of ( s j^"*^" which is useful to note for studying correlation functions is 

is Cbl + cot [^)b.p) =0. (30) 

The expression (|^ allows the reduction of the calculation of any correlation function to the calculation 
of linear combinations of quantities (s \ bp-^ . . . bp^^ \ ^' ) where the coefficients are now functions of time 
(I 5') is the initial state, and belongs to the even sector). There is an advantage in working with a quan- 
tity like (s \ bp-^ . . . bp2fj | ^'), namely its decomposition in terms of contractions of pairs, valid for some 
sets of initial states (see appendices A and B). This is particularly useful for the calculation of higher or- 
der correlation functions in those states. So, for the calculation of multi-time correlations functions like 
(s I (ii) . . . (tjv) I ^), one may use nj{t) = a|(t)aj(t) and express a](t) and a,j{t) in terms of 6],(i) 
and bp{t) using (|2j) and (p^. Then one expresses these operators in terms of 6], and bp with the help of 
( p4| ) and (|2^). Finally, one anticommutes all the fej^'s to the right and uses (pO[). Hence, one obtains a sum 
of strings which just involve the bp operators. 

The strategy for the calculation of higher order correlators is discussed in the appendices. In what follows 
we study only the time evolution of the average density at a given site j which, using the results presented 
above, is given by 

(%-W> = I E e^^"(^-^')-(-''+v)* cot (^) (. I V&-pl (31) 
p,p' 
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or, with the help of (JSC 

(f^M = 7 E e'^^'^^-^'^-^^-^+V)* (s |6t I (32) 



L 

p,p' 



4 The average density for states with random initial conditions 

Here we study the time evolution of the average density for uncorrelated random initial conditions, i.e. 
an initial state | p) where on each lattice site the probability of finding a particle is p. Such a state is a 
superposition of states with both even and odd number of particles. So, one projects it over the even and 
odd sectors, obtaining 



\eveniodd) 2 1 zb L . , /nn\ 



where ^ are the projectors over the even and odd subspace and ^ _^ (i-^2p)^ ^ normalization factor. In 
the fermion representation the state | pY^'^^ is given by 



n(l+M'cot(^)fet_^6t)|o) (34) 



I \even 2 (1 — p)^ TT /l , 2 ^ ( ''^P \ t\ 

"^^ = i + (i-2p)^ n(i + M 

where p, = p/(l - p) [|4|. 

The form ( |34|) is particularly convenient for the calculation of average values like ( s | 6p' 6p | p^"'^'^ because 
of its simple form in momentum space. If we use representation ( p9|) for ( s l"^"*^" we obtain: 

2(1 -p)^ 



> I V^pl p) 



1 + (1 - 2p)^ 



n (1 + cot ( ^ ) 6,,6_,J V n (1 + ) ^-p^^p^)l 0)-(35) 



P1>0 P2>0 



It is easily seen that if p 7^ —p' this is equal to 
2(1-P)' 



1 + (1 - 2p)i , 
Pi#p,p' 

X (0 1(1 + cot (^^^ Sp,S_,0V(l + M'cot(^^^feLp,6],,)|0) 
X (O|(l + cot(^) 6p6_p)6p(l + M'cot(^)&L^6t)|0) 
= 0. (36) 

On the other hand if p = —p', we obtain 
2(1-P)^ 



1 + (1 - 2p)L 

pi#p 



n(i..w(^)) 



X ( I (1 + cot ( ^ ) V b^P')bp' b-p, (1 + fi^cot(^) blp, bl,)\0] 



cot(^) 



1 + /i2 C0t2 ( ^ ) ■ 

So we have 



(37) 



If we insert this into (|3^) and take the thermodynamic hmit L oo, we get after a shght rearrangement 
of terms 



/)2 

-(ep+e-p)t 



'^^^^ = i/-/"^ p2 + (l-p)2tan2(|) 



(39) 



where = 1 — cos p — irj sin p. 

Notice that due to the translational invariance of the initial state there is no dependence on the site label 
j and, since ep appears in combination with e_p there is no dependence on the anisotropy parameter 77 either 
(cp + e_p = 2 — 2cosp). 

Now we keep p finite but take t ^ 00. Then, making the substitution p' = p \/t, we obtain by expanding 
the trigonometric functions 

1 p + oo -1 

<-■■"» = wfl. ^ iTS 

which is the known result given in reference Q. 

But also another scaling limit is of interest: We take t — s- 00 and p — > such that t — const. Then 
we have 

1 r~^°° ,2 1 

(hjit)) = ^ / dp'e'P -. (41) 

This integral can be calculated by the convolution theorem, giving 

{hj{t)) = pe'^''''*ET{c{2pVt) (42) 

a result also given in Notice that when py/t 00 we recover {nj{t)) oc t^^^'^. The thermodynamic 
limit ( |39| ) also holds for the odd subsector, so that for the superposition of the two states, from the even and 
the odd sector, we obtain the same result. 
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5 The local average density for the step function state 



In this section we calculate the time evolution of the local average density for an initial state where the lattice 
is empty to the left of site and full from there onwards, i.e. the initial local density is a step function. The 



fermion representation of this state, which we will represent by | <&), is, according to ( Pq ) 

L/2-1 

1$) = n 4io). 

k=0 



(43) 



The lattice sites have been renumbered from — L/2 to L/2 — 1, so that the step is located at 0, as stated 
above. Since this state is a product state in position space, it is convenient to express ( s | &J 6p' | $ ) in terms 
of ( s \alai\ $ ) with the aid of eq. (|l]) - (|22|). One gets 

7 E e^(^''-^'=)(.|ataH<I>). 



>|6t6p,|ci>; 



(44) 



k,l=-L/2 



In the infinite volume limit, L oo, ( p2| ) becomes, after the redefinition of (21), ( p^ ) such that the \/Z 
factor is absorbed into bpi , 

I ~ ' 



(45) 



Also, notice that we have absorbed the factors ^ in the redefinition of p and p', that are now continuous 
variables. Using (^), we express ( s | Sp' | $ ), in terms of its Fourier transform (s | a; | $). One can 
then express this operators back into the spin language with the help of (16h - (|lq), obtaining 



k "k+l 



\s\alai\<^>) = \ (s|s-s+|$ 



if fc = Z 



(46) 



(s |s+af+i...a^.„isj $) iffc>Z 
In this language, this expression is also easily computable for a half step profile, i.e for an initial state in 
which the lattice is empty up to site zero and half full to the right of (although the decomposition obtained 
in appendix B is not valid, because this state does not have a simple representation in terms of fermion 
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operators). The results are 



(-1)' if fc < and / > 



otherwise 



if fc < / 

if fc = ; and / > 

if fc > ; 



(47) 



for the step function profile, and 

1/4 if / > 1 and fc = 1 

1/2 / = Oandfc<0 

otherwise 



{s\alai I $1/2) = < 



1/2 



1/4 fc>land/ = fc-l 

otherwise 
for the half step profile. 

Inserting these two expressions in (^), one obtains 

00 00 00 



if fc < z 



if A; = I and / > 



if fc > ; 



(48) 



1=0 



1=0 



k=l 



for the step function profile, and 

(s $1/2) = 



(49) 



(50) 



1=0 k=l 

for the half step profile. One has to insert appropriate imaginary parts in p and p' to insure the convergence 
of the geometric series. 

It is easy to see that these expressions give the correct boundary conditions at t = 0: if one inserts ( ^9| ) 
into (^5|), the first term of the r.h.s. just gives the step function when the integrations are performed over p 
and p' . The second term is equal to X^z^o i~^y Sfe^i ^-j,k: which is identically zero. A similar result 
holds for (0). 

At finite time, one obtains from (^5|) and (^) or (|5^) respectively, the exact time evolution of the density 
profile in terms of infinite sums of products of two modified Bessel functions. More interesting, however, 
is the long time scaling limit, t 00. One performs a small momentum expansion of the trigonometric 
functions. The two expressions become identical, and we get 



1 



dp dp e 



Jpk 



(51) 



. (=0 



k=l 
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where we have performed the substitutions p — )■ p\/t and p' — > p'Vi, and accordingly, y = , I — and 
fc = The integrals can then be performed, giving 

{OO OO I 

/ — — 1 y 

We can convert the sums into integrals if we absorb a factor of t^^/^. The final result is 



{h,{t)) = i 1 + Erf(y) 



1 - Erf 



(53) 



where Erf(a:) — dte~* . A plot of ( |5^ ) is given in Figure 1. Notice that the density for large times 

depends on rj only through the scaling variable y. This means that the effect of the driving field can be 
completely absorbed by making a Galilean transformation to a frame moving with velocity rj. This agrees 
with the general result obtained in |TJ]. Also, one can easily see that this expression has the right behaviour 
when y — > ±oo. Indeed, when y +oo it agrees with the limit (^o|). When y — cx) it goes to zero as it 
should be. The diffusion front widens as but, as opposed to pure diffusion, the value of the density at 

2/ = is not the turning point of the profile function with half of the value at y = c». Instead, the turning 
point lies at a distance proportional to \/t in the initially empty region y < 0. 

This behaviour can be understood qualitatively through the following mean field argument, where no{t) 
is the density at site and noo{t) the density very far to the right of zero. Take, for simplicity 77 = and 
consider early times such that the density to the left of zero is very small and to the right is very close to 1. 
Then the mean field equation of motion for hq simplifies to: no(t) = — ^no(i) — ^nQ{t) (the second term 
of the r. h. s. is absent in the pure diffusive case). On the other hand the mean field equation for noo{t) 
is noo{t) = —n^{t). In the pure diffusive case the r.h.s. is of course zero. So the rate of decay of noo{0) 
is comparable to the rate of decay of uqIO), as opposed to the pure diffusive case in which this rate is zero. 
So, one should expect that, at large times, no{t) > ^nao{t). That is in fact what happens. Also, notice that 
the result obtained in the asymptotic large time limit is independent of whether we consider a step function 
profile or a half step function one. So we expect this result to hold for any step function with non-zero initial 
density to the right of site zero. 

Following the same mean field arguments, we should expect that if the rate of annihilation A is larger 
than 1 (A = 1 being the free fermion condition) then one would see a (transient) bump developing a,t y — 0, 
which would result from the fact that the rate for an event to the right of zero would be larger than the 
rate of an event at site at early times. It would be interesting to verify this conjecture in Monte Carlo 
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simulations. 



6 Application to experiments 

As suggested in |l^ we would like to compare our exact theoretical results for the annihilation process to 
the experimental results mentioned in the introduction by modeling the excitons dynamics by the coagulation 
process on the lattice. We have to make use of a similarity transformation to relate the two processes. 
Following @ we define a new time evolution operator by 

H = BHB^^ (54) 

where B is assumed to have the form 

B = B®^ (55) 

and the matrix B acts on a single site. Since we are studying two state systems B will be a 2 x 2 matrix. 

The diffusion-annihilation process which we are considering involves only nearest neighbour interactions. 
Hence H can be written as 

H = Y, (56) 

i 

as can be easily seen from equation (^3|) . We now investigate the relation between (^) and the Hamiltonian 
describing the diffusion-coagulation process 

AID — > IDA Dr (57) 

%A — > A% Dl (58) 

A A — y (bA (59) 

A A — > Alb 7L (60) 
with (in general anisotropic) diffusion rates Dr^l and coagulation rates ^r^l- 
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Again H is of the form (p6|). A particularly convenient way to represent and Hi,i+i is in a matrix 

form involving the two site states (there are 2x2 = 4 states, hence Hi^i+i, Hij+i are 4x4 matrices) 



i+i 



H. 



\ 






n 

u 




1-7, 

2 


2 





l-n 
2 


l+J) 
2 











A 










Dl 


-Dr 




Dl 


Dr 










7i? 



(61) 



(62) 



Since H is of the form ( p6| ) we only need to consider the effect oi B^B on (|6l]). It can be shown jl 
\ 



that the choice 
/ 

B = 



1 -1 
2 



(63) 



maps ( |6l| ) into the form ([64). The new parameters Dr, Dl, "Jr, 7l are given in terms of the old parameters 
A, 1] by 



IR = ^{>^ + V) 



(64) 
(65) 
(66) 
(67) 



Since in the exciton dynamics the rate of decay is very small compared to all other rates we can accept 
equations ( ^7| ) - (|60| ) as a good description of it. Also, as argued in ||l^, in the real system Dr = Dl = 
IR. = 7l- From equations ( |6^ ) - ( |67| ) it then follows that 77 = 0, A = 1. The first condition is just symmetric 
diffusion, but the second is the free fermion condition. Hence, our previous discussion can be used to obtain 
quantities which are of experimental relevance. However, we still have to relate the expectation values in 
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the two problems. According to ( |5^ ) the average value of some projection operator F, {s \ Fe ^* | ^ ) in 
the original problem is related to the following quantity in the new problem 



's\FB-^Be-"'B-^B\ ^] 



= (s\Fe-"'\^] 



where F = FB-^, H = B H B-^ and | ^) = S | In particular, if = hk and | = Ilz 



is a product state with average density pi on site i, then 



(68) 

/ \ 

I- Pi 



\ 



I 



F = -n. 



n 



V 



-2p, 



(69) 



(70) 



We see from (|7Cj) that for | 5') to be a physical state we must have pi < 1/2 in the original problem. 

Let us apply ( |68| ) - ( [70| ) to the case when we have random initial conditions in the coagulation problem 
with an initial density p (this corresponds to a density equal to ^ in the annihilation problem) . Then we get 
from expression ( |39| ) the following result for {fij{t))coag 



J \<^J/coag 



dp e 



-2(l-cosp)t 



p2 + (2-p)2tan2(|) 



(71) 



The scaling limit p^O, i^oo, p^t = const ( [l2|) translates to 



J \''J/coag 



pe 



P''* Eric {pVt 



(72) 



This was experimentally verified within 5 % ^ . 

The expression (fz^ ) is already known from the literature |9| where it was obtained using a continuum 
approximation for the master equation. Our calculation provides a rigorous proof of the validity of this 
expression in the scaling regime. 

A similar mapping from the step function state of section || is inappropriate because this state maps 
into a non-physical state of the coagulation problem. However, the half step function for the annihilation 
problem maps directly into the results for the step function profile in the coagulation problem and hence the 



appropriate version of (53), namely 
1 



(n,(t)). 



1 + Erf(y) 



(73) 
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can be verified experimentally. Since, as argued in the preceding section, the initial density does not seem to 
play a role, the fact that it is derived for an initial density pright = 1 is immaterial. ( |73| ) should be correct 
for any step in the late time regime. 

7 Conclusions 

We have computed time averages for some specific initial states in a simple reaction-diffusion system where 
exclusion particles hop with rates {lzL'i])/2 to the right and left respectively if the nearest sites are empty, and 
which annihilate in pairs with rate one if the nearest-neighbour site is occupied. We obtained the following 
results. 

(i) For states with random initial conditions we generalized earlier results ||^, 0, |^ to all times and 
arbitrary densities. Also, we have obtained a scaling limit for large times and small initial densities. Further, 
we have shown that due to the fact that those states are product states in the fermion representation, a 
particular product decomposition holds for quantities which are useful for the the calculation of multi-time 
and higher order equal-time correlation functions. 

(ii) For the step function and half step function state, we calculated the long time behaviour of the 
time-dependent average density. We concluded that the competition between the annihilation and diffusion 
process shifts the profile to the left of the moving front. It would be interesting to check this prediction 
experimentally. A product decomposition holds for the step function state. 

(iii) The use of the similarity transformation has shown that the scaling limit obtained above of our exact 
lattice calculation reproduces the known results obtained in a continuum approach. 

It would be very interesting to see what would happen to the density profile in the step function states if 
one goes away from the value A = 1. The results of appendix B allow for a perturbative calculation but to 
higher order in A — 1 the procedure becomes rather tedious. This could be a task that could be performed 
by simulations. 
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Appendix A Calculation of multi-time correlation functions in states 
with random initial conditions 

In section ^ we have indicated how the calculation of a multi-time correlation function could be reduced to 
the calculation of quantities like ( s | . . . 6p2„ I ^' )• Here, we show how they can be calculated for a | py^*^" 
state. 

Firstly, we define 

Z[Vp,r]-p\ = (0 |e^p>o^'°*^'^^^P^-P+'^P^P+''-^'^'P V)"""" (74) 

where the quantities rjp, rj^p are Grassmann variables obeying {r]p,r]pi} ~ 0, {rip, bp/} = 0. Their presence 
makes all the terms in the exponential commute with each other. Having this in mind and the fact that 
eVpbp = I + r]pbp {{f]p hpY = 0), we can write (|7|) as 

Z[r,p,^^p] = (0|eS.>o^°'^*'^^^-'' + Vphp)\pr^\ (75) 

p 

Notice that Z[0] = (0 | e^p>o''°*''^^^f | p)"'"" = (s|p)"^"" = 1. Also, from (|7|) and from the fact 
that rjp and bpi anticommute one gets 



dZ[rip,-q_p] 
dVpi 



- 2 

where the product over p still contains pi since bp = 0. 



One can go on differentiating with respect to the ?7p's (an even number of times). Setting rjp ^ one gets 

= (s |Sp, ...6p,„| p)''''™. (77) 

Hence Z[r]p, r]^p] is the generating function for the quantities {s \ bp-^ . . . bp^^ \ p)^^"^" . As it can be seen from 
the preceding steps, ( [fzl ) does not depend on | p)'^"'^". It is true for any state in the even sector. But the 
form (|3^ ) for | p^^"^" allows us to calculate Z[rip,ri-p] explicitly. In fact 

Z[r,p,r,.p] = AA(0|eS.>o(^°*(-^^P^-P+^P^P+'^-''^-p) J] (1 + ^2 ^ot ( ^ ) fol^, 6^, ) | ) 

p'>0 ^ ^ 

= -'^ ( I n (1 + Vpbp + V^pb^p + {cot (^^^ + rj-prip )bp b-p) 

p>0 

xn(i + M^ cot(^)bi/bi,)\o} 

= ^ J] [l + p2cot(^) (cot(^) +7?_p?7p)] (78) 



p>0 
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where J\f is the normahzation factor of | p^""^"- and the last step of ( |78| ) is similar to the manipulations done 
in equations (|3|) - (|3|). Also from Z[0] = 1 it follows that TV = [Up>a (1 + cot^i X ) ) ]"^- ^ we define 
W = In Z[r/p,77_p], we get for W[rip,ri-p] 

W[rjp,v-.p] = ^{ln[l + A*'cot(^) (cot(^)+77_p77p)] - ln[l + Ai'cot2(^)]} 

p>0 

^ / cot(f^) \ 

cot(f^) 
= 1+^.2 eot^(Y) ^-^^^ 

since {-q^pTipY = (t1-pVpY = • • • = 0. 

It is known from field theory that VF[77p,77_p] is the generating function for the connected part of the 
quantities {s\hp^ . . . bp^,^ \ pY^"^ . Since W is quadratic all the connected parts vanish for > 1 and Wick's 
theorem holds. From one concludes 

dW 

dr]p 

Also 



k,=o = 0. (80) 



d^W , /i2 cot ( ^ ) 



dvp' dvp "'^=° " 1 + cot2 ( ^ ) ^^^^ 
which agrees with (p8|). 

Further, from Wick's theorem one has 

- ( s I bp., I p)™""" ( s I 6p, I p)''"™ 

+ (^I^P.6pJp)^"'^"(5|6p,6p3|p)^'^^". (82) 

And so on. Note that as already remarked above, ( [77|) only holds for an even number of bp operators. But 
all the correlation functions can be reduced to sums of quantities like ( s | 6pj . . . bp^j^ \ p^^™ as discussed in 
section ^. 

The equality ( |8^ ) can be used to calculate ( s \ hj.^ (ti) hj^ (^2) | pY^'^'^ ■ But the manipulations necessary, 
although straightforward, become rather involved. A simplification takes place for p = 1 and p = 1/2, for in 
these cases ( s \ bpibp\ p) becomes: 

20 



For (0 = 1 

{s\bpibp\ p) = (5p,-p'; (83) 

For p = 1/2 

- - 1 / 2np' \ 



Appendix B Higher order correlation functions in the step function 
state 

In this appendix we develop a procedure for calculating higher order correlation functions in the step function 
state I $ ) defined in section |^. The quantities which we wish to compute are ( s | . . . cij^^ \ ■ Given 
that I $ ) is a product state these quantities are easily computed by expressing the annihilation operators in 
terms of the original spin operators ( ]l7| ) - ( p^ ) if ji < . . . < j2N- However, the quantities necessary for the 
calculation of correlation functions are ( s | . . . hpr^^^ \ Expressing these in terms of ( s | Oj^ . . . Oj^^ | $ ) 
by Fourier transformation involves all possible orderings of ji, . . . ,j2N and the reordering of these terms 
becomes rather cumbersome. The procedure that we develop avoids that. We start by introducing the 
coherent state 

I T^) = e-E.^fe4 I 0) (85) 

where the rjkS are a set of Grassmann variables that anticommute with the cij, aj. Developing (^5|) in powers 
of the rjk 's gives 



= I 0) - 5] 77,, 4 I 0) + 4 4 I 0) 

jl jl <j2 

jl<j2<---<jL/2 

+ •■•+ 77L/2-i---7?-L/2a^L/2---ai/2-il 0>- (86) 



The state | $ ) can be generated by taking derivatives with respect to ?7l/2-i, ■ ■ ■ iVo and then setting all 
the other 77's to zero, i.e. 



1$) = [9^0 ■•■^'ji,/2-i I '7>] all other ?7's = = {/ dr]o ■ ■ ■ dijL/2-i\ v) 



• (87) 

all other ry's = 
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The last step follows from the fact that there is no distinction between integration and differentiation in 
Grassmann algebras. Now, in order to calculate ( s | . . . aj^j^ | $ ) we define an enlarged | C) state 

l$,C) = 1/ dvo---diiL/2-ie'^^ ''^Uv)] (88) 
U J aU other ry's = 

where Q is another set of Grassmann variables anticommuting both with rjj and aj , . Using aj \ i]) — r/j \ rj) 
gives 

( I a,, . . . a,,„ I $, C > = 11 d^o . . . dr;i/2_i "^^^ ,7,, . . . r?,,„ r, ) |^ 

- 9^-,, • ■ ■ 11 d% . . . d?7L/2-i e- "^^"^ ry ) |^ (89) 

where the subscript is a reminder that we have to take the remaining 77's equal to zero. Setting the ('s 
to zero gives the desired quantity {s \aj^ . . . cij^^ \ $ ) (remember that ( s°'^'^ \ cij^ . . . cij^^ | $ ) is zero since 
I $ ) e even sector). Notice that we have not made any assumption on the order of the j's in equation (|89|), 
which may be any, as already stressed above. From (p9), we conclude that 



(90) 



is therefore a generating function for ( s | . . . Oj^^ | $ )■ 

Due to the form of (3^="^" |^ the quantity (s^"*^"! rj) can be easily calculated. We have 



(s— = (0| n(l + cot(^) 6pL,)|ry) 

p>0 

= [](l + cot(^) ,?p,y_,) = eS.x'™*^'^)''^''- (91) 
since hp\ri) = rip \ rj), where rjp is the Fourier transform of rjk ■ Substituting this in (po| ) we get 
Z[Q] ^ y d^o...dvL/2-ie-^^'^^^e^^^'°<'^^'''-^Y 



where Fki is the matrix 



sm 



^pik - I) 



(93) 



The last step of ( p2| ) follows by expressing r]p, r/^p in terms of 77^, 77; and using the Grassmann character 
of the ry's. Setting to zero all r/j's in which j < gives 
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This expression only involves the matrix elements of F such that k,l > 0, and since {s \ dj-^ . . . aj^j^ \ $ ) 
involves the derivatives of Z with respect to each ji i = 1, . . . , 2N, it will be zero if any of the j^'s is less 
than zero. 

Now equation ( p4[ ) is easily computed by completing the square in the exponential and then shifting the 
integration variables. This gives 

Z[Q] = / [ dry ] (95) 



where Fj^i is the matrix obtained by truncating Fi^i keeping only the elements k,l>0 and i^^^^ its inverse. 
The integral is standard giving 

Z[Q] = VdetFe^^"'' (96) 

It is clear that VdetF — 1 since Z\Q] = VdetF = (s| <&) = 1. To obtain the explicit form of Ff^^ we 
note that Fki is given by 

Fki = 0{k-l) - 9(1 -k) (97) 

as can be seen by Fourier transforming it, obtaining (|93|). For an L-site system ( p7| ) is an L x L antisymmetric 
matrix with every matrix element equal to minus or plus 1 depending on whether it is above or below the 
diagonal. The matrix F^i corresponds to the lower right quadrant of F^i (obtained by suppressing the first 
half of the rows and columns of F). Hence, apart from having different dimensions, F and F have the same 
simple form. This can be used to confirm explicitly that dcti^ = 1. Also, from the known Fourier transform 
of -F, or by direct methods, the inverse matrix F^^ can be shown to be 

P-^ = (-1)''-^ [9{k - I) - 9{l-k)] k,l>0. (98) 

So, we finally obtain from equation ( |9^ ) 

Z[Q] = e^^*"-' ^>^'^''^'. (99) 

As in the previous appendix this immediately shows that a form of Wick's theorem holds for the quantities 
( s I . . . aj2„ I ^' ) • For example 

( s I aj^ %3 %4 I * ) = ( s I I * ) ( s I ^3 ^4 I * ) 

- ( s I a-j, I $ ) ( s I flj, flj, I $ ) 

+ (s|a,,ajj$)(s|a,,a,3|$) (100) 
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with ( s I flji aj^ I $) = Fr^]^e{ji) e{j2) which also follows from (|99|). 
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Figure Captions 

Figure 1. Local average density profile at large times. In the horizontal axis we display y and in the 
vertical axis the density is given in units of ■ 
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